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ABSTRACT 

We report the discovery of primeval large-scale structures (LSSs) including two proto-clusters in a 
forming phase at z = 5.7. We carried out extensive deep narrow-band imaging in the 1 deg 2 sky of 
the Subaru/XMM-Newton Deep Field, and obtained a cosmic map of 515 Lya emitters (LAEs) in a 
volume with a transverse dimension of 180 Mpc x 180 Mpc and a depth of ~ 40 Mpc in comoving 
units. This cosmic map shows filamentary LSSs, including clusters and surrounding 10-40 Mpc scale 
voids, similar to the present-day LSSs. Our spectroscopic follow-up observations identify overdense 
regions in which two dense clumps of LAEs with a sphere of 1-Mpc diameter in physical units are 
included. These clumps show about 130 times higher star-formation rate density, mainly due to a 
large overdensity, ~ 80, of LAEs. These clumps would be clusters in a formation phase involving a 
burst of galaxy formation. 

Subject headings: large-scale structure of universe — galaxies: clusters: general — galaxies: formation 



1. INTRODUCTION 

Recently, there has been a great progress in observa- 
tional studies of str uctures of the high-redshift Universe. 
iSteidel et alJ l)1998f ) found a proto-cluster at z = 3.1 with 
a large population of Lyman break galaxies (LBGs). A 
significant exces s of Ly a emitters (LAEs) was found by 
iVenemans et all l)2002t) around a radio galaxy at z = 4.1, 
and a subsequent study concluded that this excess is in- 
deed a proto-cluster wi th a size of 3 — 5 Mp c |Milev et all 
I2004|) . More recently, IVenemans et alJ l|2004ft reported 
clustering of 6 LAEs near a radio galaxy at z ~ 5.2. 
For structures beyond cluster scales, wide-field narrow- 
band surveys revealed very inhomogeneous distribu- 
tion of LAEs (lOuchi et al.ll2003t IShimasaku et alJ|2003: 
lAiiki et al.ll2003t iShimasaku et al.ll2004HHu et alJl2004ji 
at z ~ 5 — 6. Shimasaku et alJ 1)200 3?) report an elongated 
overdense region of LAEs at z = 4.86 on the sky of 20 
Mpc in w idth and 50 Mpc i n le ngth. However, on the 
same sky, Shimasa ku et alJ l)2004|f found no large-scale 
structure at z = 4.79 which is closer to us by 39 Mpc, 
indicating a large cosmic variance over their surveyed vol- 
umes. Although a few segments of high-z structures have 
been found to date, the whole picture of high-z Universe 
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is veiled. Moreover, we do not know when and how the 
structures of galaxies seen at present formed from the ini- 
tial matter density fluctuations. We started a systematic 
survey to map the high-z Universe with LAEs and LBGs 
at 3 < z < 7 in the 1 deg 2 sky of the Subaru/XMM- 
Newton Deep Fie ld (SXDF: R.A = 2H 8"W, decl.= 
-5°00'00"[J2000]; ISekiguchi et"aTll200l . We obtained 
cosmic maps on scales of much larger than present-day 
LSSs to cover a network of LSSs at several redshifts be- 
tween z = 3 and 7. In this letter, we report the initial re- 
sults of our survey about structures made of LAEs at z = 
5.7. Throughout this paper, we adopt Hq = 70/i7okm s _1 
Mpc" 1 and [n m ,n A ,n,a 8 ] = [0.3,0.7,1.0,0.9]. 

2. PHOTOMETRIC SAMPLE 

We carried out extensive deep narrow-band imag- 
ing on 2003 Septemb er 27-30 and Octob er 22 with 
Subaru/Suprime-Cam l)Mivazaki et al.ll20 02L We used 
the narrow-band filter NB816 whose central wave- 
length and FWHM ar e 8150A and 12 0A, respectivel y 
HHavashino et al .120031 [Aiiki et a,l.l2003tlH 1 i et a,ll2004|K 
This filter enables us to identify LAEs, galaxies with 
a strong Lya emission line, at z ~ 5.7 ± 0.05 with a 
small fraction of foreground contaminants (20-30%; e.g. 
IHu et aLH 2004L Our 20- hour exposure imaging covered 
an area of 1.04 deg 2 with 5 Suprime-Cam pointings. The 
sky was clear and the observational conditions were sta- 
ble, with typical seeing sizes of 0". 5 — 0". 6. Th e data were 
reduced with a software package SDFRED l)Yagi et alJ 
20021 IQuchi eTafll2004al) . The sky noise of the reduced 
image is 26.0 ± 0.1 (AB) at the 5ct level with a 1".8- 
diameter circular aperture. 

We detected 3 05,012 objects down to JYB816 = 26.0 
with SExtractor ijBertin fc Arnoutal 19961) and measured 
their narrow-band excess color (i' — NB816) and contin- 
uum color (R—i') by combining exi sting very deep R- an d 
i'-band images of the SXDF (see iKodama et al.l [2004). 
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Figure^shows the two color diagram of i' — NB816 and 
R — i' for the iV£?816 detected objects, together with 
colors of model galaxies and Galactic stars. In Figure 
model LAEs at z = 5.7 have very red colors both in 
i' — NB816 and R — i' because of a strong Lya emission 
line and a continuum trough at < 1216A, which are dis- 
tinguished from colors of foreground objects. We defined 
the criteria for z ~ 5.7 ± 0.05 LAEs as i' - NB816 > 1.0 
and R — i' > 0.7, and selected 515 LAEs from all the de- 
tected objects. These selection criteria isolate emission 
line objects with an observed-frame equivalent width of 
EW ohs > 226A and a line flux of / > 6.1 x 10" 18 erg 
s _1 cm~ 2 , if a flat continuum spectrum is assumed. The 
contamination rate of our LAE sample is estimated to be 
about 30% based on our spectroscopic follow-up obser- 
vations described in section 3.2. The surface density and 
number density of these 515 LAEs are estimated to be 
£ = 0.14±0.01 arcmin" 2 and n = 5.5±0.2 x 10" 4 Mpc" 3 . 
We also calculate t he surface densit y down to the same 
magnitude limit as IHu et al . ( 2004) (NB816 < 25.05) 
a nd obtain 0.03 a rcmin~ 2 which is consistent with that 
of IHu et alJ ll200l . 

3. RESULTS AND DISCUSSION 
3.1. Primeval Large-Scale Structures at z = 5.7 

Figure His the cosmic map showing the distribution of 
515 LAEs at z = 5.7 in the SXDF. The surveyed volume 
has a transverse dimension of 180 Mpc x 180 Mpc and 
a depth of ~ 40 Mpc in comoving units. This is the first 
cosmic map, ever obtained, covering a > 100 Mpc square 
area of the Universe at any high redshifts (z > 2). In Fig- 
urc[21 the LAEs have a very clumpy distribution, forming 
concentrations with a typical size of a fe w Mpc, compara- 
ble in s ize to the proto-clu ster found by iVenemans et al.l 
l|2002|) : lMilev et all l|2004|) . These concentrations are not 
isolated, but connected with one another by filamentary 
overdense regi ons. The elongated ove rdense region of 
LAEs found bv lShimasaku et all (|2003[) at z = 4.86 may 
be a segment of a descendant of such a filamentary struc- 
ture. There are also found several voids of ellipsoidal 
shapes with sizes of 10-40 Mpc in which almost no galaxy 
exists. The characteristic sizes of filaments and voids 
seen in our map are comparable to those of the present- 
day Universe. Thus, this map marks the discovery of 
primeval LSSs at z — 5.7. 

We quantify the large-scale dumpiness of the galaxy 
distribution by estimating 020, the rms fluctuation of 
galaxy overdensity within a sphere of 20 Mpc (comov- 
ing) radius. Considering that the radial depth of the 
surveyed volume is about 40 Mpc, we estimate the fluc- 
tuation by <t 2 ~ cr| 20 = [((E 20 -S 2 o) 2 ) - S 20 ]/Si , 
where crs2o is the rms surface overdensity within circles 
of 20 Mpc radius, and X20 and £20 are the obser ved num- 
ber a nd the mean number of LAEs in a circle l|Peebles I 
1980). We obtain cr 20 = 0.4 ± 0.2. This er 2 o is compa- 
rable to or at least a half the value for the present-day 
LSSs, i.e. <72o(z = 0) = 0.5 — 0.6 ; obta ined from galaxies 
with bj < 17.5 l|Seaborne et al.lfl999|) . The characteris- 
tic shapes and the rms fluctuations indicate that these 
primeval LSSs at z ~ 6 are similar to the present-day 
LSSs. We will present results of more detailed analyses 
for these primeval LSSs, such as counts-in-cell, in our 
forthcoming paper. 



3.2. Two Clumps Identified by Spectroscopy: 
Forming Clusters? 

We find that among the dense concentrations seen in 
Figure the one at 2 h 17 m 47.2 s , -5°28'40"[J2000] has 
the highest density contrast, 6% — 3.3, with the 4.8cr 
significance level, where Ss is the surface overdensity 
for a circle of 8 Mpc radius. We refer to this concen- 
tration as Region A. Region B, a neighboring concen- 
tration at 2 h 18 m 19.6 s ,-5°32'52"[J2000], also has a high 
density contrast of <5s = 1.5 (2.2c). In order to ob- 
tain three-dimensional distributions of LAEs in Regions 
A and B, we ca rried out spectroscopy for LAEs with 
Subaru/FOCAS lIKashikawa et all2 002) on 2003 Decem- 
ber 20, 23, and 25. We placed a slit mask on each re- 
gion. Another mask was put on a control field well sep- 
arated from Regions A and B. We made a 2-hour ex- 
posure for each mask, and obtained 22 spectra in total, 
whose spectral range and resolution are A=4900-9400A 
and A/AA ~ 1000, respectively. Among the 22 objects, 
we classify one as an [Oin] emitter at z = 0.6, two as 
[On] emitters at z = 1.2, and eight as securely identi- 
fied LAEs with a characteristic asymmetric feature in 
Lya line whose skewness exceeds 1.5 with the 3<r level. 
We show spectra of these emitters in Figure |3] Since 
we have eight secure LAEs among the 11 (=8+2+1) 
identified objects, we estimate the contamination rate 
of our LAE sample to be about 30% , which is consistent 
with those of other LAE samples l|Rhoads et all 120031: 
iShimasaku et all 12003 IHu et alJ I2004D . The other 11 
objects out of the 22 are emitters with an unresolved- 
single line at ~ 8150A. Since most unresolved-single line 
emitters have been found to be L AEs with a mod erate ve- 
locity dispersion (~ 200km s-^ IHu et al1l200l . we also 
regard these 11 objects as LAEs. Thus we have 19 LAEs 
with spectroscopic redshifts as summarized in Tabled 

We plot in Figure 01 the three-dimensional distribution 
of LAEs for Regions A and B. The histogram of Figure 
0] shows that these regions have two clumps made of 10 
LAEs with a significant excess. Each of the clumps has 
a diameter of about 1 Mpc in physical units (7 Mpc in 
comoving units). For these clumps, the 1-dimensional 
velocity dispersion of galaxies is fairly small, ~ 180 
km s _1 (Clump A) and ~ 150 km s _1 (B). Although 
these clumps may not be collapsed, the formal virial 
masses can be calculated from velocity dispersion, o~ v i, 
and radius, r, with Sa^r/G, where G is the gravita- 
tional constant. We obtain ~ 1 x 10 13 M Q (Clump A) 
and ~ 8 x 1O 12 M (B). The three-dimensional density 
contrast of LAEs, S n = Sn/n, for the average of these 
clumps is ~ 80, which is comparable to those of present- 
day clusters, 100 — 200. These clumps have the high- 
est density in our surveyed volume of 9 x 10 5 Mpc 3 in 
comoving units. In the present-day Universe, this vol- 
ume typically conta ins two massive clusters wi th mass 
of 1 - 3 x 10 14 M Q l|Reiprich fc Bohringerl 12002]) . Thus, 
the discovered clumps are likely proto-clusters which are 
ancestors of today's such clusters. 

A particularly interesting feature is that these clumps 
are very high concentrations of star-forming galaxies 
with a star- formation rate of 1 — 20M Q yr _1 . The 
star-formation rate density (SFRD) in these clumps in- 
ferred from Lya fluxes is about 130 times higher than 
in the mean of the whole 1 deg 2 field, mainly due to 
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a large overdensity, 8 n ~ 80, of LAEs. This very high 
SFRD implies that these clumps would be just produc- 
ing a number of galaxies in a short period by a burst of 
galaxy formation. In contrast, present-day massive clus- 
ters are dominated by old early-type gal axies, and show 
a deficit of young; star- forming galaxies i|Dressler| 1198(1 
iButcher fc Oemlerlll98^ . Since the formation epochs of 
early- type galaxies are estimated to be z = 2 — 5 or earlier 
(e.g. IKodama. Arimoto. Barger. fc Arag'on-Salamancal 
1998), galaxies residing in these clumps are likely progen- 
itors of the old early types seen in the core of present-day 
clusters. Therefore, we are probably witnessing forming 
clusters where a number of present-day early-type galax- 
ies are just emerging as star-forming galaxies by a burst 
of galaxy formation. 

3.3. Implications for Galaxy Formation 

The Cold Dark Matter (CDM) model, the standard hi- 
erarchical scenario, predicts that very small initial fluc- 
tuations of mass density grow up gradually with time 
to evolve into galaxies and LSSs. However, the observed 
distribution of LAEs appears much more inhomogeneous 
than the dark-matter distribution, since the amplitude 
of the matter fluctuations at z — 5.7 is predicted to be 
only 1/5 of the current value. If one wants to accept 
the CDM model, a solution to this discrepancy is to as- 
sume that LAEs are biased tracers of mass in a way to 
enhance the mass fluctuations. Previous observational 
studies of high-z galax ies discuss the clusterin g bias 
iAdelberger et all 119981: iGiavalisco fc Dickinson! I200H 
IQuchi et a.U I200ll ISbimasakn et al.l 120041 IQuchi et all 
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l2004b|) . Theoretically, such biasing can be produced, if 
galaxies are formed at rare, highes t peaks of yet sma ll 
fluctuations of matter density (e.g. iMo fc Whitel [2002). 
From (T20 , we estimate that the bias parameter, 6, the 
ratio of galaxy overdensity to matter overdensity, is b = 
3.4 ± 1.8 on scales of LSSs (20 Mpc), which is compara- 
ble to the prediction at z ~ 5 from si mulations of galaxy 
formation based on the CDM model llBaugh et al.lll999t 
iKauffmann. Colberg. Diaferio. fc White! 119991) . Thus, 
the CDM model may be consistent with the early for- 
mation of LSSs that we discovered, though the model 
needs to reproduce the characteristic shapes of observed 
LSSs, long filaments and large voids. On the other hand, 
we estimate the b value, in the CDM, which reproduces 
the observed number of forming clusters (two in 9 x 10 5 
Mpc 3 ) with the size and the overdensity, 6 n ~ 8 0, fol- 
lowing the prescription of IShimasaku et aP (|2003). We 
obtain b ~ 30 for the two forming clusters, which is about 
ten times higher than the value for LSS scales, although 
this estimation of b is based on small statistics and might 
include a large error due to cosmic variance. This variety 
of b could be due to the dependence of b on the clustering 
scale or galaxy density l|Benson et alJl200H) . However, it 
is not yet clear how such a large difference in b may be 
produced within the framework of the CDM model. Our 
cosmic map may raise a question about how and where 
galaxies are formed in the early Universe. 
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Fig. 1. — Two color diagrams for continuum color (R — i 1 ) and narrow-band excess color (i' — 7VB816). The black dots indicate colors of 
305,012 objects detected with NB816 < 26.0. If the magnitude of an object is fainter than the la magnitude (la sky fluctuation), then the 
magnitude is replaced with the la magnitude. All the curves show colors of model galaxies at various redshifts. Red lines indicate model 
LAEs which are a composite spectrum of a 0.03 Gyr single burst model galaxy (GISSEL00 Bruzual & Chariot 2003) and a Lyct emission 
{EWn = 22 A); from the left to right, two different amplitudes of inter-galactic medium absorption are applied: 0.5T o ff and T e g, where r c g is 
Madau's ( 1995) median opacity. The narrow-band excess in each of the peaks in the red lines indicates the Lya emission of LAEs at z = 5.7. 
Green lines show 6 templates of nearby starburst galaxies I Kin ney et alJll996l) up to z = 2, which are 6 classes of starburst galaxies with 
E(B — V) = 0.0 — 0.7. The narrow-band excess peaks in the green lines correspond to the emission lin es of Hq (z = 0.2), [Ollll(z = 0.6) , 
H/3(z = 0.7), or [On](z = 1.2). Blue lines show colors of typical elliptical, spiral, and irregular galaxies (Coleman. Wu. & Wccdman 1980) 
which are redshifted from z = to z = 3. Yellow star marks show 175 Galactic stars given by IGunnfc Strvkcr (1983). The pink box 
surrounding the upper right region is the selection criteria of our z = 5.7 LAEs. Because the distribution of the dots show a single track 
whose colors are consistent with that of the Galactic stars, our photometry is homogeneous over the 1 deg 2 sky of the SXDF. (The paper 
with high resolution images can be downloaded from http://www-int.stsci.edu/~ouchi/work/astroph/sxds_LSS/) 
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Fig. 2. — The distribution of z = 5.7 ± 0.05 LAEs in the SXDF. The positions of LAEs are shown with yellow dots. The 
red lines correspond to contours of galaxy overdensity from <5j; = —0.25 to 3.25 with a step of A = 0.50. The characters, A and 
B, denote the positions of the two dense regions. The scale on the map is marked in both degrees and (comoving) megaparsccs. 
The large and small s quar es in the bottom l eft corner show the sizes of the surveyed areas in the previous proto-cluster searches by 
Vcncmans ct al. 12002) and Milev ct al. (2004), respectively. (The paper with high resolution images can be downloaded from 
http: / / www-int.stsci.edu/~ouchi/work/astroph/sxds_LSS/ | 
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Fig. 3. — Spectra of eight secure LAEs, together with an unresolved-single line emitter and [On] and [Om] emitters. The caption on 
each panel indicates the classification and redshift. Ticks on the vertical axis denote 3 X 10 — 19 ergs s —1 cm~ 2 A -1 , except for the following 
two panels because of display purpose. The second-top left and the fifth-top right panels have ticks of 6 X 10 -19 and 1.5 X 10~ 18 ergs s 
cm A -1 , respectively. The relative intensity of the night-sky emission is shown on the bottom right. 
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Fig. 4. — Three dimensional map of Regions A and B. The upper left panel shows the distribution of LAEs in transverse (East to 
the West) vs. radial (redshift) directions, while the bottom panel represents the distribution of LAEs projected on the sky. The red 
and blue circles correspond to the LAEs associated with the forming clusters (Clumps A and B), respectively, while the orange and cyan 
circles denote other LAEs in these regions. The white open circles plotted on the bottom panel denote the positions of LAEs without 
spectroscopic redshift. The vertical and horizontal bars indicate the length of 1 Mpc in physical units. The upper right panel shows 
the redshift distribution of LAEs with spectroscopic redshifts. The red and blue histograms correspond to LAEs in Regions A and B, 
respectively. The dotted line indicates the mean-expected number of LAEs for each region. (The paper with high resolution images 
can be downloaded from http://www-int.stsci.edu/~ouchi/work/astroph/sxds_LSS/) 
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TABLE 1 

Results of Spectroscopy 
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2.9 


02:20:21.52 


-04:53:15.0 


26.7 


5.669 


1.8 


5.1 


02:18:22.61 


-05:33:37.8 


27.6 


5.648 


1.5 


4.4 


02:18:17.35 


-05:32:22.8 


26.6 


5.642 a 


2.9 


8.2 


02:17:43.31 


-05:31:35.0 


26.5 


5.624 a 


2.5 


7.1 


02:17:45.88 


-05:27:14.5 


25.6 


5.600 


0.9 


2.4 


02:17:51.14 


-05:29:35.3 


27.2 


1.188 b 


1.7 




02:20:12.15 


-04:49:50.7 


27.4 


1.179 b 


1.5 




02:20:12.84 


-04:50:22.5 


24.5 


0.623 c 


11.9 




Note. — 
ergs s~ 1 cm" 


The quantities «^ B , / and SFR 
2 , and the star-formation rate in 


arc the i -band magnitude measured in a 2 diameter aperture, the total line flux 
units of Mq yr -1 , respectively. The redshift, z, is measured from spectra summed 


in units of 10 — 17 
over 0".9 around 



the source center, while the line flux, /, is obtained from total flux falling in the slit. 
a Securely identified LAEs 
b [Oll] emitters 
c [Om] emitter 



